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Abstract 

Monte Carlo simulations and finite-size scaling analysis have been carried out to study 
the critical behavior in a submonolayer lattice-gas of interacting monomers adsorbed on 
one-dimensional channels arranged in a triangular cross- sectional structure. The model 
mimics a nanoporous environment, where each nanotube or unit cell is represented by a 
one-dimensional array. Two kinds of lateral interaction energies have been considered: 1) 
wl, interaction energy between nearest-neighbor particles adsorbed along a single channel 
and 2) wt, interaction energy between particles adsorbed across nearest-neighbor chan- 
nels. We focus on the case of repulsive transversal interactions (wt > 0), where a rich 
variety of structural orderings are observed in the adlayer, depending on the value of the 
parameters ksT/wT (being k B the Boltzmann constant) and wl/wt- For wl/wt = 0, 
successive planes are uncorrelated, the system is equivalent to the triangular lattice and 
the well-known (a/3 x a/3) (\/3 x y/3)* ordered phase is found at low temperatures and 
a coverage, 9, of 1/3 [2/3]. In the more general case (wl/wt ^ 0), a competition between 
interactions along a single channel and a transverse coupling between sites in neigh- 
boring channels allows to evolve to a three-dimensional adsorbed layer. Consequently, 
the (\/3 x \/3) and (y/3 x y/3)* structures "propagate" along the channels and new or- 
dered phases appear in the adlayer. Each ordered phase is separated from the disordered 
state by a continuous order- disorder phase transition occurring at a critical temperature, 
T c , which presents an interesting dependence with wl/wt- The Monte Carlo technique 
was combined with the recently reported Free Energy Minimization Criterion Approach 



(FEMCA), [F. Roma et al., Phys. Rev. B, 68, 205407, (2003)], to predict the critical 
temperatures of the order-disorder transformation. The excellent qualitative agreement 
between simulated data and FEMCA results allow us to interpret the physical meaning 
of the mechanisms underlying the observed transitions. 
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I. INTRODUCTION 



Lattice-gas models have been extensively investigated in the last decades because they 
provide a theoretical framework for the description of many physical, chemical and biolog- 
ical systems. The adsorption thermodynamics and the understanding of surface phenom- 
ena have been greatly benefited from the development of these models^* 2 * 3 -. In this sense, 
the more recognized examples are the Langmuir adsorption model 2 -'- and the Ising model 
of magnetis m 2 i 5 i 6 i 7 i 8 . More recently, a number of contributions have been devoted to the 
study of adsorption of gases on solid surfac e 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 . These papers 
have included the effects of lateral interaction o 9 ! 10 ! 11 ! 12 ' 1 ? , surface heterogeneity 1 -, multisite 
occupatio n 11 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 , etc. Among them, A. J. Phares et al.— . have studied the 
structural orderings occuring in a wide variety of experimental and theoretical systems 
and its influence on the corresponding phase diagrams. 

Recently, the advent of modern techniques for building single and multiwalled car- 
bon nanotube o 23 ! 24 ! 25 ! 26 ! has considerably encouraged the investigation of the gas-solid 
interaction (adsorption and transport of simple and polyatomic adsorbates) in such a 
low-dimensional confining adsorption potentials. The design of carbon tubules, as well as 
of synthetic zeolites and aluminophosphates such as AIPO4 — 5 (Refj2&) having narrow 
channels, literally provides a way to the experimental realization of quasi one- dimensional 
adsorbents. Many studies on conductivity, electronic structure, mechanical strength, etc. 
of carbon nanotubes are being currently carried out. However the amount of theoreti- 
cal and experimental work done on the interaction and thermodynamics of simple gases 
adsorbed in nanotubes is still limite d 29130 ' 31 ! 32 ! 33 . These papers predict a complex be- 
havior. At low temperatures, Radhakrishnan et al.— have reported, by using molecular 
simulations, a phase transition occurring between two different density states, one corre- 
sponding to a low density ("gas) and the other to a higher density ("liquid) phase. The 
authors propose that this transition is induced by the interactions with the molecules 
that are in neighboring channels. A different type of transition was found by Boutin et 
alM and Lachet et al.^ at intermediate temperatures (much higher temperature than 
the phase transition considered by Radhakrishnan et al.). The molecular simulations of 
Refs . 30 ! 31 show stepped isotherms. The step was only observed using a model that includes 
both two-body and three-body lateral interactions, whereas a model with only two-body 



interactions did not show such a step. Later, this transition has been confirmed by exper- 
imental results by Martin et al.— . The step was attributed to local rearrangement of the 
adsorbed phase. Other interesting features of the carbon nanotubes have been recently 
studied. Namely, Urban et al.- 3 ^ studied the properties of an Ar film adsorbed on the 
external surface of a bundle of carbon nanotubes; and Calbi and Riccardo 3 ^ investigated 
the presence of adsorption sites and energy barriers near the ends of carbon nanotube 
bundles to determine their consequences on gas adsorption in the interstitial channels 
between the tubes. 

For theoretical purposes, adsorption in a narrowest nanotube can be treated in the one- 
dimensional lattice-gas approac h 12 ' 1 '? . This is, of course, an approximation to the state 
of real adsorbata in nanotubes, which is justified because thermodynamics and transport 
coefficient can be analytically resolved in these conditions. Many studies have been per- 
formed of specific geometries and specific adsorbate-substrate combinations^. In recent 
paper o 37 ' 38 , Trasca et al. presented simple model of a nanoporous environment. That 
is a lattice-gas model with two kinds of sites. One was a one- dimensional line of sites, 
which the authors called "axial sites, surrounded by a set of "cylindrical shell sites. The 
system was evaluated by using mean-field theoretical approach 3 - 7 - and Monte Carlo (MC) 
simulations^ 8 -. This work represents an effort in that direction. Here, we study a simplified 
lattice-gas model, which can help us to establish criteria to characterize more complex 
experimental systems. In this model, each nanotube has been represented by a one- 
dimensional chain. These chains were arranged in a triangular structure. We included 
longitudinal interactions between nearest-neighbor particles adsorbed along a single chan- 
nel, wl, and transversal energy between particles adsorbed across nearest-neighbor chan- 
nels, wt- The phase behavior depends on the values of these various energies, especially on 
the attractive or repulsive character of the interaction. In previous work, low temperature 
calculations of configurational entropy of the adlayer allowed us to identify a wide variety 
of structural orderings- 39 . Later, the influence of such structural orderings on interesting 
properties as adsorption isotherm and heat of adsorption was analyzed^. The present 
article goes a step further, studying the critical behavior of the system via MC simulation, 
finite-size scaling analysis and the recently reported Free Energy Minimization Criterion 
Approach (FEMCA) 2 ^. For this purpose, the critical temperature T c characterizing the 
transition from the disordered state to the ordered phase is obtained as a function of the 



ratio wl/wt- 

The paper is organized as follows: In Section II we describe the lattice-gas model, the 
simulation scheme and we present the behavior of T c (wl/wt), obtained by using MC 
method. In Section III we present the theoretical approach (FEMCA) and compare the 
MC results with the theoretical calculations. Finally, the general conclusions are given in 
Section IV. 

II. LATTICE-GAS MODEL AND MONTE CARLO SIMULATION SCHEME 

A. The model 

We consider the adsorption of monomers on a simple model of a nanoporous environ- 
ment. In this model, each nanotube or unit cell has been represented by a one-dimensional 
(ID) line of L adsorptive sites, with periodical boundary conditions. These chains were 
arranged in a triangular structure of size Rx R and periodical boundary conditions. Un- 
der these condition all lattice sites are equivalent hence border effects will not enter our 
derivation. The energies involved in the adsorption process are three: 

1) . Eq, interaction energy between a particle and a lattice site. 

2) . Wl, interaction energy between adjacent occupied axial sites. 

3) . Wt, interaction energy between particles adsorbed on nearest-neighbor transverse 

sites. 

Thus, the resulting substrate was an anisotropic three-dimensional array of M = L x 
R x R adsorption sites, where each site was surrounded by two "axial" sites along the 
nanotube's axis and six "transverse" sites belonging to nearest-neighbor unit cells (see 
Fig. 1 in RefM). 

In order to describe the system of N molecules adsorbed on M sites at a given tempera- 
ture T, let us introduce the occupation variable Cij± which can take the following values: 
Ci,j,k = if the corresponding site k) is empty and Cij^ = 1 if the site is occupied by 
an adatom. Then, the Hamiltonian of the system is given by, 

M 

H = w L Ci,j,kCi',f,k' + w T Cij,kCvj<,k' + (eo - y) Y Ci,j,k (1) 

(i,j,k;i',j',k') L {i,j,k;i',j',k') T i,j,k 



where (i,j,k;i',j',k') L ((i,j,k;i',j',k') T ) represents pairs of NN axial (transverse) sites 
and \i is the chemical potential. 

B. Monte Carlo simulations 

The lattice-gas was generated fulfilling the following conditions: 

• The sites were arranged in a structure of size M = L x R x R, with conventional 
periodic boundary conditions. 

• Due to the surface was assumed to be homogeneous, the interaction energy between 
the adparticles and the atoms of the substrate, eo, was neglected for sake of simplicity. 

• Repulsive transversal lateral interactions were considered, where a rich variety of 
ordered phases are observed in the planes. 

• Repulsive and attractive longitudinal lateral interactions were used. 

• Appropriate values of L and R were used in such a way that the adlayer structures 
at critical regime are not perturbed. 

In order to study the critical behavior of the system, we have used an efficient exchange 
MC or simulated tempering metho d 41 ! 42 and finite-size scaling analysi o 43 ! 44 ^ 45 . As in Ref.— , 
we build a compound system which consists of m non-interacting replicas of the system 
concerned. The m-th replica is associated with the temperature T m [or f3 m = 1/ (fcsT m ), 
being fee the Boltzmann constant]. In other words, each replica is in contact with its own 
heat bath having different temperature. Under these conditions, the algorithm to carry 
out the simulation process is the following: 

1) The compound system of m replicas is generated. For this purpose, each replica is 
simulated simultaneously and independently as canonical ensemble for ri\ MC steps 
by using a standard importance sampling MC metho d 45 i 46 i 47 . In order to determine 
the set of temperatures, {T m } ({/3 m }), we set the highest temperature, T max {(3 min ), in 
the high temperature phase where relaxation (correlation) time is expected to be very 
short and there exists only one minimum in the free energy space. On the other hand, 
the lowest temperature, T min (P max ), is somewhere in the low temperature phase 
whose properties we are interested in. Finally, the difference between two consecutive 
temperatures is set as (T max — T min ) /(m — 1) (equally spaced temperatures). 



2) Interchange vacancy-particle. The procedure is as follows: 

2.1) One of the m replicas is randomly selected. 

2.2) An occupied site and an empty site, both belonging to the replica chosen in 2.1), 
are randomly selected and their positions are established. 

2.3) By using a standard Kawasaki algorithm^, an attempt is made to interchange 
the occupancy state of the sites chosen in step 2.2). 

3) Exchange of two configurations X m and X m i, corresponding to the m-th and m'-th 
replicas is tried and accepted with the probability W (X m , /3 m \X m i, /3 m /). In general, 
the probability of exchanging configurations of the m-th and m'-th replicas is given 

\ 1 for A < 

W(X m ,/3 m \X ml ,/3 ml ) = \ (2) 

I exp(-A) for A > 

where A = (/3 m — /3 m /) [H(X m i) — H(X m )]. As in Ref.— , we restrict the replica- 
exchange to the case m m + 1 . 

4) Repeat from step 2) m x M times. This is the elementary step in the simulation 
process or Monte Carlo step (MCS). 

The procedure l)-4) is repeated for all lattice's sizes. For each size, the equilibrium 
state can be well reproduced after discarding the first MCS. Then, averages are taken 
over n MC s successive MCS. The canonical expectation value of a physical quantity A is 
obtained in the usual way as follows: 

i nucs 

( A K = - — E M* m (t)] (3) 

timcs t =i 

All calculations were carried out using the parallel cluster BACO of Universidad Na- 
tional de San Luis, Argentina. This facility consists of 60 PCs each with a 3.0 GHz 
Pentium-4 processor. 

As it is standard for order-disorder phase transitions, a related order parameter was de- 
fined. In particular, at 9 = 1/3 [2/3] (being 9 = N/M the surface coverage ), a(\/3x\/3) 
[(\/3 x v^3)*] ordered structure is formed in the planes below the critical temperature. 
Depend on the sign of the longitudinal interactions, the order is propagated to all planes. 



For repulsive wl, adatoms avoiding configurations with nearest-neighbor interactions or- 
der along the channels in a structure of alternating particles separated by empty sites. On 
the other hand, attractive monomer-monomer longitudinal interactions favor the forma- 
tion of pairs of nearest-neighbor adsorbed particles along the nanotubes. The resulting 
structures are shown in Figs. 5 and 8 of Ref.-^. 

For the case of repulsive longitudinal interactions, Fig. 1 a) shows two successive planes, 
k and k + 1, for one possible configuration of the phase appearing at critical regime and 
9 = 1/3. Due to the periodic boundary conditions the degeneracy of this "local phase" 
is equal to six. These configurations allow us to decompose the "local lattice" into six 
different sublattices [see Figs. 1 b)-c)]^. The coverage on each sublattice is denoted as 
9 s (s = 1, 6). In this way, an "local order parameter", ipk, can be defined as 

Vk= E %-o t \ (4) 

where we sum the differences (in absolute value) between the coverage corresponding to 
two sublattices. 

When the system is disordered (T > T c ), all sublattices are equivalents and the order 
parameter is minimum. However, when a configuration of the local phase appears at 
low temperature (T < T c ), this is allocated on a sublattice. Let us suppose that this 
configuration lies on the sublattice s. Then, the coverage 9 S is maximum (9 S = 1) and 
the coverage of the rest of the sublattices is zero or minimum. Consequently, <pk is also 
maximum. 

On the basis of <fk, the generalized order parameter, ip, can be written as, 

L 

v = ^E^ ( 5 ) 

k=0 

where A is a normalization factor. The definition (jSJ) is computationally convenient and 
ip appears as a good order parameter evidencing the order-disorder phase transition. 

In a similar way, it is possible to define the order parameter corresponding to 9 = 2/3 
and repulsive longitudinal interactions. 

For attractive longitudinal interactions, it is not appropriate to define local sublattices. 
In this case, each sublattice s(s = 1, 3) lies on the total lattice (see Fig. 2) and ip can 
be easily defined as ip = \9i — 9 2 \ + \9 1 — 9 3 \ + \9 2 — 9 3 \ 



Now, the reduced fourth-order cumulant, Ur[l], introduced by Binder— and related to 
the order parameter, can be calculated as: 



where Ur [Ul] represents the cumulant obtained by variable R [L] and fixed L [R]. The 
thermal average (...)r, in all the quantities, means the time average throughout the MC 
simulation. 

The standard theory of finite-size scalin g 43 i 44 i 45 allows for various efficient routes to 
estimate T c from MC data. One of these methods, which will be used here, is from the 
temperature dependence of Umn(T), which is independent of the system size for T = T c . 
In other words, T c is found from the intersection of the curve Umn (T) for different values 
of R[L], since Ur[l](T c ) = constant. 

C. Computational results 

The thermodynamic properties of the present model have been investigated by means 
of the computational scheme described in the previous section. As a consequence of the 
equivalence particle- vacancy, the critical behavior at 9 = 2/3 is as at 9 = 1/3. Then, we 
restrict our calculations to 9 = 1/3. In addition, we set wt = 1 and vary wl/wt from — 1 
to 1 (-1 < w l /wt < 1). 

In order to understand the basic phenomenology, we consider in the first place null lon- 
gitudinal interactions {wl/wt = 0). In this particular case, successive planes are uncorre- 
lated and the system is equivalent to the well-known triangular lattice. The value obtained 
of kBT c /wT = 0.3354(1) confirms this arguments and validates the MC schem o 50 i 51 i 52 i 53 . 
The data are not shown here for brevity. 

Hereafter, we discuss the behavior of the critical temperature as a function of wl/wt- 
We start with the case of attractive longitudinal interactions. As an example, Fig. 3 
illustrates the reduced four-order cumulants plotted versus k B T/wx for wl/wt = — 1. 
From their intersections one gets the estimation of the critical temperature. The lattice 
sizes used in the simulation^ 4 - are compiled in Table I along with the values of the param- 
eters in the simulated tempering runs. In the figure, the critical temperature is obtained 
from the curves of Ur{T) (calculated for different values of R and fixed L). The resulting 



value, ksTc/wT = 0.7817(1), agrees very well with previous determinations reported in 
the literature^. In Ref.— , a value /tbT c /% ~ 0.76 was obtained from the inflection on the 
function s(T), being s(T) the configurational entropy per site of the adlayer as a function 
of the temperature. Due to the finite-size scaling technique and the simulation procedure 
used in this contribution, the estimation of T c in the present work is expected to be more 
accurate than that reported previously. 

The study was extended to repulsive longitudinal interactions. Fig. 4 shows the data 
for a typical case (wl/wt = 1), resulting k B T c /wT = 0.6098(5), in well agreement with 
the value obtained in Ref.^ [k B T c jwT ~ 0.59). As indicated in Fig. 3, the parameters of 
the simulation are listed in Table I. 

Due to the presence of anisotropy (the couplings are taken to be different in the dif- 
ferent lattice directions), it is expected that the correlation functions in transversal and 
longitudinal directions may be governed by correlation lengths diverging with different 
critical exponent a 55 ! 56 . However, it is worth pointing out that we do not assume any 
particular value of the critical exponents for the transitions analyzed here in order to 
calculate their critical temperatures, since the analysis rely on the order parameter cu- 
mulant's properties^.. In addition, the procedure shown in Figs. 4 and 5 was repeated 
for the curves of Ul{T), which were obtained for variable L and fixed R. As an example, 
Fig. 6 presents the results obtained for the case wl/wt = — 1. As it is expected, identical 
results (within numerical errors) are obtained in both ways. 

Finally, the calculations were carried out for wl/wt = 
-0.75,-0.50,-0.25,0.00,0.25,0.50 and 0.75 and the results were collected in Table 
H*&. As it can be observed, the critical temperature presents an non-trivial behavior 
as a function of wl/wt- An understanding of the dependence of fcsT c /w T on wl/wt 
can be developed by following the subtle interdependence of energetic and entropic cost 
necessary to alter the ordered phase. This will be discussed in Section III. 

III. THEORETICAL APPROACH: FREE ENERGY MINIMIZATION CRITERION 

Hereafter, we will use FEMCA2S in order to discuss the dependence of ksTc/wT vs. 
wl/wt obtained from MC simulation. 

In a closed system of adsorbed particles with repulsive interactions, the phase transition 
occurring in the adsorbate is a continuous (second-order) phase transition. In other words, 



the entropy, S, varies continuously from a completely ordered state (when T —>■ 0) to a 
disordered state (when T — > oo). Around T c , S changes abruptly (but continuously)^. 
Then, it's possible to analyze the phase transition taking into account the Helmholtz free 
energy, F = E — TS (being E the mean energy), in the two extreme states (maximum 
order and maximum disorder). Accordingly, 



= lim F and F = lim F (7) 



then 



F oo «F T>T C (8) 

F oo »F => T<T C (9) 

^ = ^0 => T^T C (10) 

The last equation allows to determine T c . This calculation is not exact due to the system 
does not pass from an extreme order to an extreme disorder. There exist intermediate 
states between the two extreme states. However, as we will show in the following analysis, 
the eq. (jl()j) provides a very good approximation for T c . Interested readers are referred 
to Ref.— for a more complete description of FEMCA. 

In general, for a system of interacting particles at temperature T results: 

/o = e - Ts and = - Ts^ (11) 

where e and s represent the mean energy per site and the entropy per site in the thermo- 
dynamical limit, respectively, 



E S 
e = lim — and s = lim — (12) 

M->oo M M^oo M 



If /o = foe, this is 



e - Ts = eoo - Tsoo 



(13) 



then T ~ T c and 

m Ae — en ,„ 

T c « — = — (14) 

As Soo - s 

From eq. ()14|). it is possible to calculate the critical temperature and to interpret the 
dependence of UbTc/wt with wl/wt obtained from simulations. As in Section II, we 
restrict the study to 6 = 1/3(2/3) and —1 < wl/wt < 1. 

Case I: 9 = 1/3 and wl/wt > 

In general, e^^) can be calculated from mean- field approximation. Thus, 



eoo(#) = ^ (GN6w T + 2Ndw L ) (15) 



In this case = N/M = 1/3, and 



eoo(l/3) = ^w T + ^w L (16) 

In order to calculate the entropy of the disordered state, the configurational factor of 
monomers Q, is employed 



_ M\ 

~ N\{M-N)\ (17) 

Thus, 



Sqo = hm — — — [ll 

M^oo M 

In the particular case of 6 = 1/3, the entropy per site of the disordered state results, 



Soo (l/3) = -k B (\n^ + hn2^j (19) 

In addition, the mean energy per site and the entropy per site for the ordered state at 
6 = 1/3 and T = are eo(l/3) = so(l/3) = 0. Then, the critical temperature depends on 
the mean energy and the entropy of the disordered state. From eqs. (jTHJ), (fT^j) and (fT3|). 
we obtain T c (l/3): 



r e (i/3) 



eoo(l/3) 
,(1/3) 



+ K (2Q) 



-k B (ln| + |ln2) 



^ Tc(1/3) « K > and w L /w T > 0) (21) 



Finally, 

ks' ^ 

Wq ^ — In I — I In 2 

Case II: 9 = 1/3 and wl/wt < 



The mean energy and the entropy of the disordered system are as in eqs. (II 5|) and (|18j). 
On the other hand, for the ordered system, sq(1/3) = and eo(l/3) = wl/3. Then, 



^5^(1/3) 3 IZt 



~ , , , , , (% > and w L /w T < 0) (22) 

— In | — ~ In 2 

As it is expected, the calculations for 8 = 2/3 (do not shown here) provide identical 
results as Cases I and II. 

Fig. 6 shows the comparison between the simulated results previously presented in Table 
II and the theoretical predictions obtained from FEMCA for the critical temperature 
as a function of wl/wt- The MC simulations reveal the main characteristics for the 
behavior of the critical temperature versus Wl/wt' i) the curve presents a minimum for 
wl/wt = 0; and ii) for negative values of wl/wt, the critical temperatures are higher 
than the corresponding ones for positive wl/wt's- Both characteristics are very well 
reproduce by FEMCA. 

The physical meaning of the main features of the critical temperature can be interpreted 
from the theoretical approach. In this framework, the eq. (|T4|) shows that ksTc/wT 
depends on the mean energy and the entropy of the disordered state. The behavior of 
these quantities as a function of wl/wt allows to understand the arguments presented in 
the previous paragraph. Thus, the values of Soo, so and are identical for repulsive and 
attractive longitudinal interactions. In addition, the magnitude of eo is constant (eo = 0) 
for wl > and increases with wl (e = wl/3) for wl < 0. Now we can interpret the 
difference between the two regimes in Fig. 6. From —1 < wl/wt < 1, the variation in 



the entropy is constant. On the other hand, the variation of the mean energy increases 
linearly with Wl, being higher for Wl < than for Wl > 0. 

IV. CONCLUSIONS 

In the present work, we have addressed the critical properties of a simple lattice-gas 
model, which mimics a nanoporous environment, where each nanotube or unit cell is rep- 
resented by a one-dimensional array. The results were obtained by using MC simulations, 
finite-size scaling theory and the recently reported FEMCA, which is based on a free 
energy minimization criterion. 

The system was characterized by two parameters wl/wt and ksT/wT, being wl and wt, 
the longitudinal and transversal energy, respectively. We focused on the case of repulsive 
transversal interaction energy among adsorbed particles (wt = 1) and 6 = 1/3(2/3), in 
such a way that a rich variety of ordered phases are observed in the adlayer: 

• For wl/wt = 0, the system is equivalent to the well-known triangular lattice in 2D. 

• For wl/wt < 0, the formation of pairs of nearest-neighbor adsorbed particles along 
the nanotubes is favored. Consequently, the (\/3 x \/3) and (\/3 x a/3)* phases are 
reinforced and extend along the channels. The critical temperature decreases from 
0.7817(1) for w L /w T = -1 to 0.3354(1) for w L /w T = 0. 

• For wl/wt > 0, the (y/3 x y/3) and (y/3 x y/3)* structures are formed in the planes 
at low-temperatures and order along the channels in a array of alternating particles. 
The critical temperature increases from 0.3354(1) for w L /w T = to 0.6098(5) for 
w L /w T = 1. 

With respect to the analytical approach, FEMCA provides results in very good quali- 
tative agreement with MC simulations and constitutes a theoretical framework in order 
to interpret the behavior of ksTc/wx vs wl/wt in the critical concentrations. 

Future efforts will be directed to (a) include attractive wt longitudinal interactions 
between the adparticles, (b) obtain the phase diagram A; B T c /w T versus 9 in the whole range 
of coverage, and (c) develop an exhaustive study on critical exponents and universality. 
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V. TABLE AND FIGURE CAPTIONS 



Table I Parameters of the simulated tempering runs for two typical cases (wl/wt = —1,1)- 

Table II Critical temperatures corresponding to the critical coverage 9 = 1/3(2/3). The data 
were obtained from the crossing of the cumulants. 

Figure 1: a) Snapshot of two successive planes, k and k + 1, for a possible configuration of the 
ordered phase appearing at 9 = 1/3 and repulsive longitudinal interactions. Solid 
circles represent occupied sites. The different sublattices used in order to define a 
local order parameter characterizing this low temperature structure are shown in 
parts b)-c). 

Figure 2: Different sublattices defined for attractive longitudinal interactions and 9 = 1/3. 

Figure 3: Ur(T) and Ul(T) [inset] versus ksT/wr, for a typical case of Wl < 0: 9 = 1/3 and 
w L /w T = -1. 

Figure 4: U R (T) versus k B T/w T for 9 = 1/3 and w L /w T = 1. 
Figure 5: Ul(T) versus ksT/wr for 9 = 1/3 and wl/wt = — 1. 

Figure 6: Comparison between simulated and theoretical results for k B T c /wT vs. Wl/wt at 
9 = 1/3(2/3). The dotted lines are a guide for the eyes. 
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